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Composition modulation in Pb(Sc,Nb)0 3 alloys * 



L. Bellaiche and A.M. George 
Physics Department, University of Arkansas, 
Fayetteville, Arkansas 72701, USA 

A first-principles-derived approach is used to study properties of Pb(Sc.Nb)C>3 alloys that 
are compositionally modulated along a direction away from the total electrical polarization. 
The direction of the local electrical polarizations (centered in the planes perpendicular to 
the modulation direction) is found to be spatially-dependent, and can be adjusted with the 
modulation parameters. Furthermore, modulated structures with large composition gradi- 
ents adopt a new monoclinic phase for their ferroelectric ground state, and exhibit very 
large piezoelectric and dielectric coefficients. These unusual features can be viewed as a 
consequence of a composition modulation-induced internal electric field. 

* Supported by ASTA grant N99-B-2L NSF grant DMR-9983678 and OXR grant 
N00014-00- 1-0542. 
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EXHIBIT 




Examples of recent discoveries : 



• 1) High piezoelectricity in single-crystals 
of PZN-PT and PMN-PT (~ 2,000 pC/N) 

[S.-E. Park and T.E. Shrout, J. Appl. Phys. 82, 1804 (1997)] 

• 2) Monoclinic phase in the morphotropic 

phase boundary of PZT 

[B. Noheda et al, Appl. Phys. Lett. 74, 2059 (1999)] 

• 3) Asymmetry of the double-well potential 

in AB0 3 / AB'0 3 / AB"0 3 structures 

[N. Sai et al, Phys. Rev. Lett. 84, 5636 (2000)] 

Ferroelectric alloys show great promise for 
the development of new materials 



Aim of the present study : 

To look for a mechanism: 

1 . Occurring in ferroelectric alloys 

2. Yielding new structural properties. 

3. Leading to huge electromechanical responses 

-► This mechanism is the composition modulation 

Plan of the talk : 

L Method 

II. Composition modulation in Pb(Sc,Nb)0 3 

III. Understanding the unusual features 



I. Method: Heff for A(B Vx B" x )0 3 

- Mimic the true alloy A(B'i-x B" x )0 3 
by a virtual simple A<B>0 3 system 

— E,({ti},u) 

depends on strain {r| } and mode u 

[ Zhong et al, PRL 73, 1861 (1994) ] 

- Consider alloying as a perturbation 
of the A<B>0 3 system 

-> E = E, ({n}, u) + E 2 ({ti}, u, {a} ) 
Incorporates alloy conf. {a} 

[ Bellaiche et al, PRL 84, 5427 (2000) ] 

- All the parameters of Ei and E2 are 
determined by first-principles calculations 

H 



Expression of E2 in Pb(Sc,Nb)0 3 



O: = +1 : Niobium atom on site j 
Gj = -1 : Scandium atom on site j 

Most important interaction: Intersite effect 
between atom on site i and local mode on site i 

E 2 = Ei ( Lj Q j? jOj e „ ). Ui = - Zj (Sj e „ ). Z* u { 



all negative unit vector 
(up to 3 rd nn) from j to i 



Nb on site j site i u . 



"Electric-field* 
generated by 
atom sitting on 





Accuracy of Heff for Pb(Sc 0 5 Nb 0 5 )0 3 



1. Disordered Pb(Sc 0 . 5 Nb 0 5 )O 3 (PSN): 
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Correct phase transition and good agreement with 
experiment [* Chu et al, J. Appl. Phys. 77, 1671 (1995)] 

2. Rocksalt-ordered PSN at small temperature: 

Heff LDA 
Rhombohedral with Rhombohedral with 
<u>=0.0675(x +y +z) <u> =0.0679(x +y +z) 



II. Composition Modulation 
in Pb(Sc 1 . x(k) Nb x(k) )0 3 

The composition x(k) varies along [001] 



x(k) = 

t 

Integer indexing 
the (001) planes 



x 
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Average 
composition 
(50 %) 



A f(k, X) 



Shape 



Amplitude 

Wavelength 



Atomic disorder is assumed in each (OOl)plane 



Sinusoidal modulation with A=0.25 and X=12a 
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□ □ 



□ ux(k) = uy(k) 
^-uz(k) 

— <ux> = <uy> = <uz> 
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Plane Index k 



Rhombohedral ground state (Total P //[ 1 1 1 ]) 
uz(k) sinusoidally varies as a function of k 
► Local polarizations are spatially-dependent 



Saw modulation with A=0.25 and A,=12a 
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d ux(k) = uy(k) 
-^uz(k) 

— <ux> = <uy> = <uz> 
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Rhombohedral ground state 

The local polarizations can be shaped by 
adjusting the parameters of the modulation 



Is it possible to change the ground-state? 
Saw modulation with A,=4a and tuning A 

(SC o.5 Nb o.5 ) / (SC o. 5+ A Nb 0.5.A ) ' 

(Sc 0.5 Nb 0.5 ) / (Sc o. 5 . A Nb 

0.5+A J 



0.08 




Amplitude A 

- Rhombohedral phase for small A (P //[ 1 1 1 ]) 

- Monoclinic phase for intermediate values of A 

► P between the [1 1 1] and [110] directions 

- Orthorhombic phase for the largest A (P // [1 10] ) 
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Effect on piezoelectricity (at 50K)? 
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Huge shear piezoelectric coefficient! ! ! 

(around the monoclinic-orthorhombic transition) 



Effect on dielectric response (at 50K)? 



7000 




Amplitude A 



Large dielectric susceptibility 

(around the monoclinic-orthorhombic transition) 



III. Understanding the unusual features 

Composition modulation generates an internal 
electric-field along [001], acting on the k th plane: 

E (k) = g({ Qj J) [x(k+l)- x(k-l)] 

Strongly negative 
in PSN 

Modulation parameters 
play a role 

E(k) < 0 when x(k+l) - x(k-l) > 0 
► uz(k) < <uz> 

E(k) > 0 when x(k+l) - x(k-l) < 0 
► uz(k) > <uz> 

E(k) = 0 when x(k+l) - x(k-l) = 0 
— ► uz(k) = <uz> 



Sinusoidal modulation with A=0.25 and A,=12a 
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Effect on local modes at 50K 
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Local polarizations are shaped by 
the internal electric-field!!! 



Why is it possible to alter the ground-state? 
Saw modulation with A,=4a and tuning A 



Planes Electric-field 

(Sc 0 . 5 Nb a5 ) 

(SC 0.5+A Nb 0.5-A ) 

(Sc 0.5 Nb ) 

(SC 0.5-A Nb 0.5+A ) 

Magnitude of E is linearly dependent on A 
— ► <uz> decreases when increasing A 

— ► Monoclinic phase for intermediate A 
— ► Orthorhombic phase for the largest A 





Conclusions 

- Use a first-principles derived scheme to study 
properties of compositionally modulated PSN 

- These structures exhibit unusual features 

► Local polarizations that can be shaped 

by adjusting the modulation parameters 

— ► New structural phases 

► Large electromechanical responses 

- These unusual features result from the 
existence of an internal electric-field 
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Polarons and Oxygen Vacancy Clusters in Ferroelectric Thin Films: Vacancy 
Self -patterning and "Fatigue as a Phase Transition" 

J. F. Scott 
Cambridge University 



We have carried out infrared absorption studies of polarons in both lead 
zirconate titanate (PZT) and barium- strontium titanate (BST) as functions 
of temperature (20 - 650 K) , donor and acceptor doping, fatigue, and 
carrier concentration. The random fields associated with Ba/Sr occupancy 
in BST or Zr/Ti site occupancy in PZT facilitate Anderson localization. 
Oxygen vacancies order in [ID] chains and/or [2D] rafts. Anderson 
localization in dimensions less than [3D] have a zero-field threshold. 
Thus, strong polaron localization is expected in such ferroelectric films. 
In agreement with this we find strong polaron absorption in BST and PZT 
near 4000 cm-1 but not in pure strontium titanate or SBT (strontium bismuth 
tantalate) . The polaron mass M** is found to be 16.m(elec.) = 3 M* (band 
mass). This may affect transport and leakage current models. Second 
harmonic generation is nominally forbidden in BST and SrTi03 but occurs 
strongly. Its temperature dependence (increasing by x5 from ambient to 4 
K) gives a measure of oxygen vacancy cluster size and is fitted to the 
detailed theory of Prosandeev; this is a no- free-parameter fit, with the 
constants obtained independently from the temperature at which the 
reciprocal dielectric constant extrapolates to zero from above (36±3 K) and 
the lowest symmetric q=0 phonon energy at 4 K (ca. 40 cm-1) . The results 
are related to the concept of "fatigue as a percolative phase transition", 
and it is found that Steinsvik's oxygen vacancy activation enthalpy H in 
strontium titanate varies near a critical concentration of Xc = 7% as 

g 

H(x) = (Xc-x) 

where g = 7/4 is the percolation exponent for a cubic lattice. At Xc = 7% 
a phase transition to a 2x2x1 superlattice occurs. The oxygen vacancies 
order along [101] axes according to Becerro et al . ; this is the axis of the 
main domain walls involved in pinning and fatigue (Janovec, private 
communication) as observed by McKinstry, Randall, et al . in PZT. 



POLARIZATION DYNAMICS IN FERROELECTRIC THIN FILMS AND NANOSTRUCTURES 

R. RAMESH, Department of Materials and Nuclear Engineering and Center for 

Superconductivity Research, Department of Physics, Materials Research 

Science and Engineering Center, University of Maryland, College Park, MD 

20742. 

Over the past several years, there has been considerable focus on 
ferroelectric thin films, specifically for possible use in nonvolatile 
memories. A significant portion of the work to date has focused on thin 
film processing, process integration and reliability studies, specifically 
approaches to solving the problems of fatigue and imprint. However, in 
order for thin films of PZT and SBT (the two primary materials systems under 
consideration) to be successfully implemented in nonvolatile memories, 
several other basic and applied problems have to be resolved. Specifically, 
the switching dynamics of the spontaneous polarization under applied field; 
the long term relaxation of the remanent polarization state and the scaling 
of the ferroelectric properties with lateral and vertical dimensions are 
issues of a fundamental nature and require rigorous understanding. From 
another point of view, ferroelectric and dielectric perovskites present very 
fascinating systems where the dynamics of the spontaneous polarization or 
that of dipolar clusters over the time scale of lpsecond to 10 years 
(almost 20 decades of time) is of great relevance from both the fundamental 
and applications point of view. The use of novel surface probes such as 
electric force microscopy, scanning capacitance microscopy is now enabling 
high resolution studies of the ferroelectric and dielectric properties. In 
this presentation, I will present results of studies aimed at understanding 
the interplay between thin film processing, size effects, film 
microstructure and polarization dynamics (specifically relaxation dynamics) . 
Sub-micron ferroelectric capacitors have been fabricated by focused ion beam 
milling and are being probed using a combination of conventional 
characterization techniques as well as novel surface probes such as scanned 
force microscopy and spectroscopy. This approach has enabled us to image 
polarization relaxation phenomena at a resolution of 20nm ; these 
measurements coupled with theoretical treatments have revealed the details 
of nucleation and growth of reverse domains. We have been able to directly 
measure the ferroelectric and piezoelectric properties of sub-lOOnm 
capacitors. In such sub-micron structures we find that the polarization 
dynamics of 90=B0 domain walls is very different from continuous films. 
Specifically, we find that by releasing the clamping effects imposed by the 
substrate through focused ion beam milling, 9O=B0 domain walls can be moved by 
the application of an electric field. This observation is of critical 
importance in enabling approaches that can take advantage of the large (4%) 
strain that is generated by the movement of the 90=BO domain wall. We have 
also found that we can control the location and the polarization direction 
within the 90=B0 domain through the use of vicinally cut substrates. The 
results of these experiments point to approaches to create novel, 
electrically tunable ferroelectric gratings. 



The work at Maryland is supported by the NSF-MRSEC. 



Domain Walls and Compositional Disorder in Pb-based Perovskites 



D. Vanderbilt, 1 B. Meyer. 1 R. Hemphill. 2 L. Bellaiche, 2 and A. Garcia 

1 Department of Physics and Astronomy, Rutgers University 
2 Physics Department, University of Arkansas. Fayetteville 
Departamento de Fisica Aplicada II, Universidad del Pais Vasco, Spain 



Two aspects of recent work on the properties of Pb-based ferroelectric 
perovskites will be discussed. First, a thorough understanding of ferroelec- 
tric switching requires a microscopic description of the underlying domain 
walls and their dynamics. We have carried out ab-initio calculations of 90° 
and 180° ferroelectric domain walls in PbTi0 3 . We find that it is possible to 
construct metastable supercell structures representing both kinds of domain 
walls using supercells of only modest size. The relaxed structure, energy, 
and polarization profile are calculated for each case. Two features of the 90° 
dislocation are surprising: the barrier for its motion is found to be miniscule, 
and a small offset of the electrostatic potential across the domain wall is 
found. Second, a first-principles-derived approach is used to study the prop- 
erties of rocksalt-ordered and disordered Pb(Sc 0 .5Xb 0 .o)O 3 (PSN) alloys [1]. 
The paraelectric-to-ferroelectric transition temperature T c is strongly depen- 
dent on the atomic configuration, while the piezoelectric response vs. T/T c 
is nearly independent of the chemical order. Our calculations are consistent 
with the experimental finding [2] that ; at T = T c . ordered PSN undergoes 
a normal ferroelectric transition, while disordered PSN transforms from a 
relaxor state to a ferroelectric state. 

[1] R. Hemphill. L. Bellaiche. A. Garcia, and D. Vanderbilt, Appl. Phys. 
Lett. 77, 3642 (2000). 

[2] F. Chu ; I.M. Reaney, and N. Setter, J. Appl. Phys. 77, 1671 (1995). 



Polarization Rotational, Elastic, and Electromechanical Properties of Piezoelectrics 

Ronald E. Cohen and Eugene Heifets, Carnegie Institution of Washington, 5251 Broad Branch 
Rd., N.W. Washington D.C 20015 and California Institute of Technology, MS 252-21, Pasadena 
CA 91 125 and H. Fu, Rutgers University 

Macroscopic polarization in crystals is a vector quantity that is a bulk property, and the 
conventional picture of maximum coupling to electric fields was to apply an electric field parallel 
to the polarization to induce a strain, in order to maximum the coupling between the electric field 
and polarization vectors. However, the electromechanical response is governed by the 
piezoelectric constants, which are third-rank tensors, and in general are not required to give a 
maximum strain when the polarization and applied field are parallel. Fu and Cohen [1] proposed 
that this effect is responsible for the huge piezoelectric response found in PZN-PT and PMN-PT 
piezoelectrics [2] through polarization rotation. This apparently has become quickly accepted, 
and new experimental data supports this idea. For example, an intermediate monoclinic phase has 
been discovered in PZT between the rhombohedral and tetragonal phases [3], that is understood 
as the rotation of the polarization from (1 1 1) to (001), so apparently polarization rotation even 
occurs in the absence of an applied electric field. The idea is very simple: rhombohedral 
perovskite ferroelectrics (e.g. the ground state of BaTi0 3 have a small strain, whereas the lattice 
strain is much larger in tetragonal ferroelectric perovskites (1\% in BaTi0 3 and 6% in PbTi0 3 , for 
example). By applying a field along (001), oblique to the (1 1 1) rhombohedral polarization, one 
can obtain a huge strain as the polarization is rotated to (001). This can also be thought of as a 
field induced phase transition from rhombohedral to tetragonal, through an intermediate 
monoclinic phase. In some cases this may be a continuous transition, and in other cases there may 
be a first-order transition at some applied field, and thus hysteresis and a discontinuity in strain 
with field. 

The computation of macroscopic polarization using the Berry's phase approach [4] is the key to a 
first-principles understanding of piezoelectricity. We have used the Linearized Augmented Plane 
Wave (LAPW) method within the LDA and GGA to compute the piezoelectrics constants of 
PbTi0 3 and PZT 50/50 [5]. For PbTi0 3 we find good agreement with experiment, though there is 
a wide spread among experimental measurements. There are three pieces to understanding 
electromechanical response, piezoelectric constants, elastic constants, and dielectric constants. 
We are now studying the elastic properties of PZT, which should be useful since PZT single 
crystals have not been successfully grown, so single crystal constants are not directly accessible 
experimentally. 

This research is supported by the Office of Naval Research. Computations were performed on the 
Cray SV1 at the Carnegie Institution of Washington, supported by NSF and the Keck Foundation. 

[1] H. Fu and R.E. Cohen, Nature, 403, 281-283 (2000); H. Fu and R.E. Cohen, In: Fundamental Physics of 
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[3] B. Noheda et al., Appl. Phys. Lett. 74, 2059, 1999; B. Noheda et al., Phys. Rev. B 61, 8687, 2000. 

[4] R.D. King-Smith and D. Vanderbilt, Phys. Rev. B, 47, 1651, 1993; R. Resta, Rev. Mod. Phys., 66, 899, 

1994. 

[5] G. Saghi-Szabo, R.E. Cohen, and H. Krakauer, Phys. Rev. Lett. 80, 4321, 1998; 59, 12771, 1999. 



Lattice Distortions in Pb(Zr,Ti)0 3 Alloys Near the 
Morphotropic Phase Boundary 

David J. Singh and Marco Fornari 
Center for Computational Materials Science 
Naval Research Laboratory 
Washington, DC 20375 

First principles calculations within the framework of density functional theory have 
shown that both ferroelectric distortions that provide piezoelectric couplings and instabili- 
ties consisting of octahedral rotations are close in energy over a wide range of compositions 
in the rhombohedral part of the Pb(Zr,Ti)0 3 phase diagram. Additionally based on the 
opposite pressure dependencies of these two basic instabilities, it has been suggested that 
they may coexist locally near the morphotropic phase boundarj\ This is due to the cou- 
pling to local strain fields arising from the B-site cation disorder. Here we review these 
results and present details of first principles structural relaxations for supercells with a 
variety of local Zr-Ti orderings. These are to show this effect, as well as to elucidate the 
Pb displacements for various local chemical environments in the disordered Pb(Zr.Ti)0 3 
alloy The results are discussed in terms of the Pb-0 chemistry particularly hybridization 
of Pb states with O p states and the interplay between displacive and rotational instabil- 
ities, particularly the fact that variations in the amplitude of the ferroelectric instability 
can occur on shorter length scales that that of rotational instabilities in Pb(Zr ? Ti)0 3 . 

*Supported by the Office of Naval Research 



"Local atomic structure of PbZr03-PbTi03 near the morphotropic phase boundary". 

W. Dmowski, T. Egami , L . Farber and P. K. Davies 
Department of Materials Science and Eng., 
University of Pennsylvania 
Philadelphia, PA 19104-6272 

We have studied the local atomic structure of PbZr03-PbTi03 solid 
solution (P2T) near the morphotropic phase boundary (MPB) using pulsed 
neutron scattering. PZT is applied as solid state sensors and actuators 
that couple strain to electric signal. On the Zr-rich side PZT is 
antiferroelectric and has the orthorhombic Pbam structure up to 7% of PbTi03 , 
while above 7% it becomes ferroelectric with the rhombohedral structure. The 
PZT phase diagram is characterized by the almost vertical MPB [1] at 
48-52% of PT, which separates the Zr rich rhombohedral phase from the 
tetragonal (P4mm) Ti rich phase. At this phase boundary the piezoelectric 
response of PZT is exceptionally high and stable over a wide range 
of temperature. Recently a high resolution x-ray diffraction study of PZT 
at the MPB at low temperatures revealed a previously unknown monoclinic 
phase Cm [2] . However, it is difficult to determine local atomic dsplacements 
away from the high symmetry positions from the crystallographic studies. 
In our study, in addition to conventional crystallographic 
Rietveld analysis, we used the pair distribution function (PDF) method. 
Since the PDF is obtained by the direct Fourier- trans formation of both the 
Bragg peaks and diffuse scattering intensity the inter-atomic distances 
are obtained without the assumption of the unit cell symmetry. 
While the crystal structure changes substantially through the MPB, we found 
that the local atomic structure does not show dramatic changes. 
The actual change is more of a statistical nature in the population of 
Pb displacements. Crystallographically the Pb displacement 
is in [100] in the tetragonal phase and [111] in the rhombohedral phase, 
with the rotation from one to the other in the intermediate monoclinic phase. 
Locally, however, pure [111] Pb displacements are not observed even in the 
Zr-rich rhombohedral phase. The [100] displacements are majority in the 
tetragonal phase, while they are gradually replaced by the [110] 
displacements through the MPB. The distortions of the Ti and Zr oxygen 
octahedra are very close to those observed in pure PT and PZ. The observed 
gradual change in the local structure through the MPB suggests that the 
"polarization rotation" is driven by the gradual change 

in the partition of the Pb displacement patterns, and the crystal structure 
simply reflects the change in the average structure. 



Synchrotron X-ray Studies of Superlattice Ordering in Lead Magnesium Niobate Single 

Crystal Doped with 6.25% Ti 

A. Tkachuk, P. Zchack, E. Colla and Haydn Chen 

Department of Materials Science and Engineering 
And Frederick Seitz Materials Research Laboratory 
University of Illinois, Urbana, Illinois 61801 



Abstract 



Temperature dependence of the correlated displacements in small regions of 
approximately 30 A has been studied in the lead magnesium niobate (PMN) single crystal 
containing 6.25 % Ti (PMN-PT), using synchrotron x-ray scattering techniques. Local 
phase transition temperature has been found by monitoring superlattice reflections 
located at the face-centered locations in the reciprocal space, known as a spots. These 
superlattice reflections are related to the atomic displacements, which are correlated on 
the short-range scale in the nanoregions that are different from the chemical 
nanodomains. Displacements are believed to be anti-ferrodistortive in nature and their 
magnitude decreases with increasing temperature below the freezing temperature, Tf - 
260 K, while the radius of these nanoregions was found to be temperature independent. 
Tf=260 K, that marks the local structural phase transition in the nanoregions, is very close 
to the phenomenological freezing temperature, obtained from the Vogel-Fulcher fit to the 
frequency dependence of the dielectric response for this particular composition. A 
comparison of PMN-PT results to the pure PMN, where nanodomains due to both anti- 
ferrodistorions and chemical ordering were found, will be presented. Structure-property 
relationship in the relaxor systems will be also discussed. 

* Work supported by the US Department of Energy 



Broadband-Brillouin Scattering of Relaxor Ferroelectric Crystals 
Seiji Kojima and Fuming Jiang 

Institute of Materials Science, University of Tsukuba, Tsukuba, Ibaraki, 
305-8573, Japan 

The fast relaxation process of relaxor ferroelectric crystals is studied by 
broadband-Brillouin scattering technique. The quasielastic scattering 
appears at about 100K above the cubic to tetragonal transition temperature 
Tct in a PMN-PT crystal. On further cooling its intensity increases down 
to Tct . Since the spectral resolution is very high, the relaxation time 
is determined accurately as a function of temperature. Consequently, the 
fast relaxation in GHz region is clearly observed. Similar quasielastic 
scattering is also observed in PZN-PT and SBN crystals. The change of sound 
velocity and damping is also observed above Tct. By combining a microscope 
with a Fabry-Perot Interferometer, the position depnednece of elastic 
properties is also studied by changing the spot area to be observed. 



Improving effective Hamiltonians for simple and complex ferroelectric systems 
Karin M. Rabe, Rutgers University 

First-principles effective Hamiltonians have proved generally quite 
successful in simulating the finite-temperature structure and 
properties of ferroelectrics and related materials. However, there are 
quantitative discrepancies with experimental results, especially in 
transition temperatures. These are in part due to the use of 
approximate density functionals and to idealization of the real 
systems, and in part due to the additional approximations associated 
with the effective Hamiltonian approach. In this talk, I will discuss 
a number of issues relevant to the latter, including choice of optimal 
prototype structure, the construction of improved lattice Wannier 
functions, relaxations in solid solutions and thin films, 
"high-accuracy computation of the quadratic part of the effective 
Hamiltonian, and higher-order anharmonic terms such as the scaling of 
interatomic force constants with distance. Examples will be drawn from 
the ferroelectric semiconductor GeTe, bulk and thin film ! PbTiO$_3$, 
and PZT. 



Composition modulation in Pb(Sc,Nb)0 3 alloys * 



L. Bellaiche and A.M. George 
Physics Department, University oj Arkansas, 
Fayetteville, Arkansas 72701, USA 

A first-principles-derived approach is used to study properties of Pb(Sc,Nb)0 3 alloys that 
are compositionally modulated along a direction away from the total electrical polarization. 
The direction of the local electrical polarizations (centered in the planes perpendicular to 
the modulation direction) is found to be spatially-dependent, and can be adjusted with the 
modulation parameters. Furthermore, modulated structures with large composition gradi- 
ents adopt a new monoclinic phase for their ferroelectric ground state, and exhibit very 
large piezoelectric and dielectric coefficients. These unusual features can be viewed as a 
consequence of a composition modulation-induced internal electric field. 



* Supported by ASTA grant N99-B-2L NSF grant DMR- 9983678 and ONR grant 
N00014-00-1-0542. 
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Domain patterns and macroscopic electro-mechanical behavior 



Kaushik Bhattacharya 
California Institute of Technology 
Pasadena, CA91125 

We report on our recent theoretical and experimental work on understanding the electromechanical 
behavior of ferroelectrics, with particular attention to link between macroscopic properties and 
domain patterns. We develop a continuum theory based on energy minimization, similar to the 
classical Landau-Devonshire theory, and use it to understand the domain patterns in single crystals, 
polycrystals and thin films. We study the effect of applied stress and electric field, and use the 
theory to find an electro-mechanical loading path that yields large electrostriction. We validate 
these results experimentally by demonstrating strains close to 1% in BaTi0 3 single crystals under 
constant mechanical compression and cycling electric field. In-situ microscopy shows that strains 
are generated by domain switching mediated by domain wall propagation. We conclude with some 
recent theoretical results on polycrystals: (i) why it easier to pole PZT near the morphotropic phase 
boundary and (ii) the effect of texture on electrostriction in thin films. 



Cation ordering in [Nax/^Bii^TiO^ 
Prediction of the ground-state configuration. 

B.P. BURTON and E. COCKAYNE 
Materials Science and Engineering Laboratory, Ceramics Division 
National Institute of Standards and Technology, Gaithersburg, MD 20899 

The Vienna Ab-initio Simulation Package (VASP) was used to perform 
fully relaxed planewave pseudopotential calculations of formation energies 
(AEvasp) for a large number of ordered supercells in the perovskite based 
system NaTiO^ - BiTiOs; including 30 supercells at the 
[Nai/ 2 Bii/ 2 ]Ti03 composition. ^Prom the set, {AE VA $p}, the lowest 
energy configuration, and presumed ground state (GS) is a 40 atom 
supercell (40GS?) with space group symmetry PAifmmc (neglecting 
ferroelastic distortions). Its structure is characterized by pseudocubic 
doubling of the primitive perovskite cell, plus alternating [100] rows of Na 
and Bi atoms in (hkO) planes and alternating [010] rows in (hk^) planes: 
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One expects this structure to be easily disordered because [^00] or [0^0] 
stacking faults can be introduced without changing any cation-cation pair 
correlations for pairs closer than fourth nearest neighbors. Accordingly, it 
is no surprise that AE§%j$> = 163.17 kJ/AB0 3 , is almost identical to 
that for an 80-atom supercell, AEvasp = 163.02 kJ/AB0 3 , which is 
related to GS40? by a [^00] stacking fault, and doubling of the c-axis. 

A succession of cluster expansion (CE) Hamiltonians were fit to 
{AEvasp}, and then used to calculate the set {AEce} for all possible 
configurations of the 80-atom cell. For each CE, {AEce} contained a new 
candidate GS, so AE$^ P was calculated, added to {AEvasp}, and a 
new CE fit was performed. No candidate GS, however, had 
ae vasp < & e vaTp- Therefore GS40? is the presumed GS for 
NaTiOs - BiTiOz, but more CE convergence tests are required to confirm 
this. 

1 



Geometric Formulation of Quantum Stress Fields in Density Functional Theory- 
Christopher L. Rogers and Andrew M. Rappe, 
Department of Chemistry, University of Pennsylvania 

The quantum-mechanical stress is a fundamentally important quantity 
which provides insight into the relationship between the electronic 
structure and the mechanical properties of condensed-matter systems. 
Nielsen and Martin have provided a tractable and well-defined 
expression of the macroscopic or total stress for periodic solid-state 
systems within the framework of density functional theory (DFT) , 
allowing the study and prediction of a large class of structural 
properties from first principles. 

The quantum stress field provides significantly more information than 
the total stress, allowing one to study the spatial dependence of the 
stress and its relationship with local electronic and structural 
phenomena. In particular, the quantum stress field can provide new 
insights into the polarization and piezoelectricity of perovskites. 
However, previous derivations of the quantum stress field are not 
unique since they are determined only up to a gauge term. 

We will present a unique expression for the quantum stress field 
within the local density approximation (LDA) . Our formulation 
exploits a known relationship between the strain tensor field and the 
Riemann metric tensor, giving rise to a quantity which is gauge 
invariant with respect to choice of energy density. Stress field 
results from DFT- LDA calculations of solid molecular hydrogen will be 
reported. Extensions of the formalism to pseudopotentials and gradient 
corrected exchange-correlation functionals will be discussed. 



Dynamics of Relaxors 



Robert Blinc, R. Pirc, and V. Bobnar 
J. Stefan Institute, Ljubljana, Slovenia 



Relaxor ferroelectrics are characterized by a high value of the static dielectric 
constant and a strong frequency dispersion of the complex dielectric permitivity. In 
analogy with dipolar glasses the longest relaxation time diverges at a freezing 
temperature 7}. It is now ge nerally accepted t hat the propert ie s of relaxors are due to 
nanogiz gd polar clusters ,. We have proposed a spherical random bond - random field 
(SRBRF) model, which is capable of explaining a number of static properties of 
relaxors. The basic reorientable dipoles in relaxors are polar clusters, which vary in 
size and orientation. In view of that the psedo spin which is proportional to the dipole 
moment of the polar cluster is not a fixed length vector as in dipolar glasses, but can be 
regarded as a continuous variable subject only to the spherical constraint. Here we 
present a dynamical SRBRF model based on the interacting polar cluster picture. We 
assume that the reorientation of polar clusters is the leading relaxation mechanism 
above the freezing temperature. The Langevin equations of motion for the polar cluster 
pseudo spins are written down and an exact stationary solution is obtained in the 
asymptotic time limit from which the linear and non-linear dielectric permitivities are 
obtained. The theoretical predictions are compared with experimental results and a 
theoretical pressure-temperature phase diagram of relaxors is proposed. 



Ferroelectricity in doped incipient 
ferroelectrics. 



Yizhak Yacoby and Yakov Girshberg. 

Racah Institute of Physics, Hebrew University, Jerusalem, Israel 91904 



Abstract. 

We present a comprehensive model of ferroelectricity in incipient ferroelectrics doped 
with ferroelectric ions such as KTaC>3:Nb. We show that a strong electron - longitudinal 
phonon interaction at the impurity sites together with the presence of virtual electrons 
in the conduction band of the host lattice are responsible for the spontaneous off- 
center displacements of the impurity ions in the cubic phase. The interaction of these 
off-center displacements with the host lattice soft mode drives the ferroelectric phase 
transition. We show that taking the tunneling of the off-center ions into account we can 
predict quantitatively the critical impurity concentration for T c =0, the concentration 
dependence of T c , The temperature dependence of the dielectric constant and the 
temperature dependence of the soft mode. Finally, taking spin-spin correlations into 
account we find that the system has a central peak very close to T c and the soft mode 
saturates at a finite value of a few cm" 1 . 



What happens inside a polarized dielectric? 

Paolo Umari, 1 Andrea Dal Corso, 2 and Raffaele Resta 3 



^IRRMA-Institut Romand de Recherche Numerique en Physique des Materiaux, 

CH-1015 Lausanne, Switzerland 
2 INFM-Scuola Internazionale Superiore di Studi Avanzati, 
Via Beirut 4, 1-34014 TYieste, Italy 
3 INFM-Dip. di Fisica Teorica, Universita di Trieste, 
Strada Costiera 11, 1-34014 Trieste, Italy 

We address the very basic issue of what happens, at a microscopic level, inside a 
polarized dielectric. We show that the complete information about electronic polarization 
is embedded in the microscopic polarization Pind(r)- In a crystalline dielectric this vector 
field is defined by its Fourier coefficients Pg- Previous studies in the literature have 
addressed the induced electronic charge density (alias the divergence of our vector field) 
where the most relevant information is obliterated. In fact, the coefficients of the induced 
charge are iG ■ Pq^ while the macroscopic polarization is simply the G = 0 coefficient P 0 . 
This shows once more that macroscopic polarization has nothing to do with the periodic 
charge of the polarized dielectric: such fact goes against a common belief, shared by many 
textbooks, and based on the Clausius-Mossotti model. 

We provide a quantum-mechanical expression for P ind (r), and we present first-prin- 
ciple results for two case studies: an ionic crystal (NaBr), and a covalent one (Si). The 
vector field has the dimensions of a dipole per unit volume: we don't find, however, useful 
to visualize it as a "dipole density" . thinking in terms of a continuous distribution of 
dipoles. Instead, the physical meaning of P (ind) (r) is better understood by imagining that 
the applied field is adiabatically switched on in time: P (ind) (r) is then proportional to 
the microscopic current flowing through the sample while the field is switched on and the 
dielectric is polarized. In a Clausius-Mossotti model such transient current does not cross 
the cell boundaries, whereas in a real crystal it percolates through the sample. We show 
that our case studies are both far from the Clausius-Mossotti limit. In the case of Si, 
the (unperturbed) valence charge defines a continuous network of bonds. When a field is 
switched on, most of the polarization current Pj n d(r) flows within narrow channels along 
the bonds, and steers to keep clear of the core regions: the dominant feature is that the 
current percolates across the material. Although less dominant, a similar feature occurs 
even in NaBr: a big surprise, since , ionic crystals are usually regarded as prototypical 
Clausius-Mossotti materials. 

For ferroelectrics, which are mixed ionic/covalent materials, we expect a dominant role 
of percolating transient currents. The present viewpoint is likely to provide a thorough 
microscopic understanding of the main trends (i.e. high or low polarizability) in different 
materials. 

Partly supported by ONR. 



Lattice dynamics of PMN measured by inelastic neutron scattering. 



S. Gvasaliya, S. Lushnikov, S. Vakhrushev 
A.F. Ioffe Phys.-Tech. Institute, 194021 St. Petersburg, Russia 

B. Dorner, A. Ivanov 
Institute Laue-Langevin, F-38042 Grenoble, Cedex (, France 

D. Strauch, K.Schmalzl 
Theoretische Physik, Universitat, D-93040 Regensburg, Germany 

Despite long history of study the details of lattice dynamics of the relaxor ferroelectrics remain 
unclear. There are sill hot discussions in the literature about the origin of the Raman spectra in the 
relaxors [1]. Recently performed neutron inelastic scattering measurements revealed the soft 
excitation [2] and disappearance of the low TO mode at small q [3]. To get a "complete" picture 
of the phonon dispersion curves in the model relaxor crystal PMN inelastic neutron scattering 
measurements for energy transfer up to 25 THz were carried out at 12K (to reduce usual 
anharmonic multiphonon processes) for main high-symmetry directions and lattice-dynamical 
shell-model calculations were performed. Obtained results are much more complicated than those 
observed in simple perovskites. There are 4 branches in [100] direction with well defined 
polarization at small q, below 6.5 THz a broad unresolved continuum between 6.5 and 12 THz, 
well defined excitation at 13.4 THz observed at conditions favorable for the observations of 
longitudinal as well as transverse phonons, and the weak broad signals above 16THz. 6.5 to 12 
THz band was observed in longitudinal and transverse geometries and in many Brillouin Zones 
(BZ), but we were not able to discover any trend such as to follow a strong signal throughout one 
BZ Low energy excitations were fitted to the shell model for the virtual crystal with B cation 
having the mass equal to 1/3 m Mg +2/3m Nb . No calculated counterpart was found for the 13.4 THZ 
excitation. Some possible ideas to explain the observed results will be discussed. 



1. See e.g. I.G.Siny, S.G. Lushnikov, R.S. Katiyar, V.H. Schmidt, Ferroelectrics 226, 
191 (1999). 

2. A. Naberezhnov, S. Vakhrushev, B. Dorner, D. Strauch, H. Moudden, Eur. Phys. J. 5, 
11,13 (1999). 

3. P.M. Gehring, S.B. Vakhrushev, G.Shirane "Fundamental physics of ferroelectrics 
2000" ed. R.E. Cohen, AIP Conference Proceedings, 535, 314 (2000). 



Importance of quenched-disorder upon the structure-property relations 
of oriented piezocrystals 

Dwight Viehland 

Naval Seacommand, Newport, RI 02841 

Investigations of oriented PMN-PT crystals have been performed by various 
methods. The results have demonstrated the importance of quenched 
disorder upon the unique properties of these materials. Several important 
features have been found these are: (i) the observation of a significant 
ferroelastic switching under uniaxial stress; (ii) various intermediate 
metastable conditions; (iii) the breakdown of uniform polarization under 
continuous unipolar excitation; and (iv) the observation of polarization 
relaxational dynamics which scale with the Random-field theory under 
bipolar drive during polarization reversal. These data are explained using the 
concept of symmetry adaptive ferroelectric mesostates. Various symmetry 
adaptive states exist including monoclinic and orthorhombic, between the 
stable rhombohedral and tetragonal phases. Symmetry adaptation occurs due 
to accomodation of the large transformational strain. 



17 0 Isotropic Chemical Shift Probing of Phase Transitions 
in Hydrogen-bonded Solids: Squaric Acid 

K. L. Pierce, R. Fu and N. S. Dalai 
Department of Chemistry and National High Magnetic Field Laboratory 

Florida State University 
Tallahassee, FL 32306-4390 

O-H'-O and N-H-O type bonds are known to be involved in ferroelectric phase 
transitions and other cooperative phenomena shown by a variety of solids, but the 
underlying mechanisms is not yet well understood. Current theoretical models assume 
relatively little or no role for the O atom, but this is perhaps related to the lack of suitable 
experimental data, in particular that obtainable by high resolution NMR. Here we present 
our 17 0 MAS results on squaric acid, H2C4O4, through its paraelectric-antiferroelectric 
phase transition at 100°C. The data clearly demonstrates that MAS using single crystals 
provides much (-factor of 4) narrower peaks than with powders, pointing to the role of 
the anisotropic bulk magnetic susceptibility (ABS) broadening mechanisms [KJymachyov 
and Dalai, Z Phys. B 104, 651, 1997), Fu et al.,7. Phys. Chem. 102, 8732 (1998)]. The 
spectrum consists of four well resolved peaks at T<T N (-373K), and coalesces to a singlet 
above T>Tn- Several novel features are notable, including the coexistence of the spectra 
from the paraelectric and an ti ferroelectric phase in a short temperature regime around 
373K, two 17 0 sites in the paraelectric phase and the dynamic coalescence of the spectra 
at 373. The latter is in contrast to earlier NQR results where the coalescence was noted to 
occur at about 425 K. There is a large change in the 17 0 istotropic chemical sift at T N , 
indicating that the phase transition must involve a displacive component, and not solely 
the order-disorder behavior as assumed in general. 



First Principles study on the coexistence of magnetism and 
ferroelectricity in YMnO$_3$ 

Alessio Filippetti and Nicola A. Hill 
Materials Research Lab, University of California at Santa Barbara. 



Hexagonal YMnO$_3$ is one of the rare materials known to be both 
ferroelectric and magnetic. The occurence of these two phase transitions 
in the same compound stimulated the idea of designing new, powerful 
technological devices. However, the coexistence of magnetism and 
ferroelectricity has been object of not much theoretical study so far, 
and in this work we attempt to understand how magnetic and electric 
polarizations are related to each other at microscopic level, and which 
conditions have to be fulfilled for the coexistence to be realized. 
Specifically, we present results for YMnO$_3$ (ferroelectric and 
antiferromagnetic) , BaTiO$_3$ (ferroelectric and non magnetic), and 
CaMnO$_3$ (non ferroelectric and antiferromagnetic) . Our study reveals the 
fundamental role of the hybridization between cation d and oxygen p 
orbitals in f erroelectrics : a strong change of s-p charge hybridization is 
always associated (and should be considered the driving force) to 
ferroelectric distortions. Since changes of hybridization are large for 
materials with mostly unoccupied d orbitals, the possibility of having 
magnetic f erroelectrics seems to be ruled out. However, this impediment 
can be overcome in compounds where the electric polarization is 
constrained to a specific direction. Indeed, only orbitals oriented along 
this direction contribute substantially to the ferroelectric distortion, 
whereas orbitals orthogonal to the distortion do not need to fulfill the 
d$"0$-ness criterion, thus they can carry a net magnetic moment. According 
to our calculations, these conditions are verified in hexagonal YMnO$_3$: 
the electric polarization is oriented along the c axis, and a magnetic ' 
moment equal to 3.7 $\mu_B$ is calculated on each Mn ion, consistent with 
a S=2 spin state. Furthermore, we find that non magnetic YMnO$_3$ can not 
be ferroelectric, thus for YMnO$_3$ magnetism not only does not prevent, 
but is also necessary for the ferroelectricity to occur. 



First-principles study of Pb 2 MgTe0 6 perovskites 



Razvan Caracas and Xavier Gonze 
University Catholique de Louvain, Unite de Physico-Chimie et de Physique de Materiaux, 1, place Croix du Sud } B-1348 

Louvain-la-Neuve, Begium 

Pb2MgTe06 perovskite, elpasolite, is cubic, Fm3m, at high temperature and becomes incom- 
mensurately modulated with a rhombohedral, #3, average structure at lower temperatures. Both 
the cubic and the average structure have been investigated using the density functional theory. El- 
pasolite is an isolant with a 1.5 eV gap. The electronic band structure is done of weakly dispersive 
electronic bands. The peaks in the electronic density of states are assigned to the different atomic 
and/or molecular orbitals. The ab initio determination of the structure is in good agreement with 
the experiment. The phonons in T are analyzed in terms of group theory and compared with the 
experiment. The Born effective charges and a phonon band structure are also reported. 
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Why is there an isotope effect on T c in KH 2 P0 4 upon deuteration but not when replacing 

16 Oby 18 0? 

Annette Bussmann-Holder and Naresh Dalai 
Max-Planck-Institut fiir Festkorperforschung, Heisenbergstr. 1 , D-70569 Stuttgart, Germany 
Department of Chemistry and NHMFL, Florida State University, Tallahassee, FL32306, USA 

High-precision NMR experiments on KH2PO4 (KDP) crystals and its analogues have revealed 
that the ferroelectric phase transition observed in the KDP-family carries a pronounced 
displacive component, quite opposite to the general believe that this transition is a 
prototypical order-disorder driven transition. In consensus with former work the huge isotope 
effect on T c upon deuteration has been confirmed and it was observed that also the deuterated 
compound shows clear evidence that a coexistence of order-disorder and displacive dynamics 
are present. Most interestingly oxygen isotope experiments have not been carried through yet 
to search for an isotope effect on T c stemming from the oxygen. Especially the recent finding 
of oxygen isotope induced ferroelectricity in SrTiC>3 could have lead to speculate that a 
pronounced enhancement of T c should also be detected in KDP. Yet opposite to these 
assumptions recent new I8 0 NMR measurements on KDP revealed that an isotope effect on 
T c practically does not exist. This finding can be explained within a coupled pseudospin- 
lattice-electron model, and supports strongly our previous findings of coexistence of order- 
disorder and displacive dynamics in hydrogen-bonded systems. In addition it is shown, why 
different results are obtained in the quantum paraelectric SrTi0 3 , but should not be obsrevable 
in the analoguous compound KTaC>3. 



MODEL OF POLAR CLUSTERS IN RELAXORS : 
CHARGE TRANSFER AND LOCAL CONFIGURATION INSTABILITY 

EFFECTS 

V.S. Vikhnin 

A.F. Ioffe Physical Technical Institute, St. Petersburg 194021, Russia 
E-mail: valentin.vikhnin@pop.ioffe.rssi.ru; 
valentin.vikhnin@physik.uni-osnabrueck.de 

A model of polar clusters (PC) and of chemical fluctuation clusters (CFC) in 
PMN-like relaxors is proposed. The model is mainly based on the following 
phenomena. First is a charge transfer state formation corresponded to (i) 
localized Charge Transfer Vibronic Excitons [1-4] induced by Mg-Nb disorder 
on the one hand, and to (ii) electronic and hole polaron appearing due to local 
charge compensation in CFC on the other. Second is Local Configuration 
Instability phenomenon which is realized in relaxors both in the cases of 
electronic Jahn-Teller polaron states as well as of cluster ordering states. Local 
Configuration Instability effects lead to (iii) crossover between small and 
intermediate electronic polaron states in CFC regions, to (iv) Charge Trasnsfer 
Vibronic Exciton local ordering [1-3] in PC, and to (v) antiferrodistortive local 
transformation in CFC. The origin of co-operative phenomena, of Giant 
Piezoelectric Effect (GPE), and of characteristic dielectric losses including 
Vogel-Fulcher law validity in relaxors finds possible explanation on the basis of 
microscopic approach. For instance, GPE is induced by an appearance of mixed 
TA - soft TO lattice modes as well as of mixed CFC soft antiferrodistortive, and 
PC soft polar quasilocal modes in the vicinity of PC due to giving rise a strong 
local piezo-electric effect. Here two addition cross-type lattice-cluster mixed 
modes are also active in the GPE formation. The drastic increasing of local 
piezo-electric effect due to the contribution of the soft antiferrodistortive 
quasilocal CFC mode is induced by the proximity to the Local Configuration 
Instability point. Recent experimental results of Anomalous Synchrotron X-ray 
Scattering (A. Tkachuk, H. Chen, P. Zschack, and E. Colla, 2000), of EPR (N. 
Takesue, Y. Fujii, K. Koyama, M. Motokawa, and H. You, 2000), and of NMR 
(R. Blinc et al, 1999 - 2000) investigations can be also explained within the 
model proposed. . 

Supported by German DAAD, and by Russian RFBR Programs. 

1. V.S. Vikhnin. Ferroelectrics, 1997, v. 199, p.25. 

2. V.S. Vikhnin. Z. Phys. Chem., 1997, v.201, p.S201. 

3. V.S. Vikhnin. Ferroelectrics Letters, 1999, v.25, p.27. 

4. V.S. Vikhnin, R. Blinc, R. Pirc. Ferroelectrics, 2000, v.240, p.355. 
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Pressure as a Probe of Ferroelectric Properties: Quantum Regime* 



G. A. Samara and E. L. Venturing Sandia National Laboratories 
L. A. Boatner, Oak Ridge National Laboratory 

Quantum fluctuations can strongly influence the low temperature response of a system near a 
structural phase transition. Among the manifestations of quantum fluctuations at ferroelectric 
phase transitions are different critical exponents, the suppression of the transition temperature, 
T c , below its classical value and the ultimate development of a quantum paraelectric state. To 
study these quantum effects as well as the crossover from the classical to the quantum regime, it 
is necessary to shift T c continuously from high to low temperatures. This is conventionally done 
by chemical substitution, but high pressure is a "cleaner" variable. Results on ABO3 perovskites 
exhibiting normal ferroelectric and relaxor behavior at 1 bar will illustrate the above quantum 
effects at high pressure. Emphasis will be on KTaC>3 and KTN crystals doped with < 0.05at% Ca 
or Ba. The addition of these two dopants introduce dipolar defects into the KTa03 lattice which 
strongly enhance the dielectric susceptibility of the host lattice in the quantum paraelectric state. 
In a Ca-doped KTN crystal with 2.3at% Nb pressure induces a crossover from normal 
ferroelectric to a relaxor state which on further increase in pressure crosses over to a quantum 
paraelectric phase. The susceptibility of the paraelectric phase is very strongly dependent on 
pressure and biasing electric field. The results will be discussed in terms of the physics involved. 

♦This work was supported by the Division of Materials Sciences, Office of Basic Energy 
Sciences, United States Department of Energy under Contracts No. DE- AC04-AL85000 at 
Sandia National Laboratories and No. DE-AC05-96OR-22464 at Oak Ridge National 
Laboratory. 



Bound Charge Diffusion and Polar Nanocluster Dynamics in Proton Glass Crystals 

V. HUGO SCHMIDT, Physics Dept., Montana State Univ., Bozeman, MT 59717 
Slater proposed that KH 2 P0 4 (KDP) has 2 polar and 4 nonpolar (with higher energy E b ) 
H2PO4 H-bond configurations. Takagi postulated H3PO4 and HPO4 groups with still 
higher energy E a . Via intrabond H transfer, Takagi pairs form, diffuse independently, 
and recombine or else annihilate with new partners, allowing H-bond reconfiguration. In 
both KDP-type ferroelectrics and in proton glasses such as Rbi, x (NH 4 ) x H2P04 (RADP), at 
higher temperatures there are many Takagi groups. They have short life and are likely to 
annihilate with new partners, leading to nearly Debye-like dielectric relaxation. Below 
the ferroelectric transition in KDP, each Takagi pair member requires energy Eb to 
diffuse one step farther from its partner. This energy bias leads to their rapid 
recombination. In proton glasses with intermediate x, there is no ferroelectric or 
antiferroelectric transition. Well above the freezing range, the Takagi concentration falls 
rapidly with temperature and the Takagi group lifetime and path to annihilation grow 
long. The energy landscape of the annihilation path must be level on the averaged 
However, the random energy steps ±Eb or 0 make the highest barrier encountered be 
proportional (on average) to the square root of the diffusion path length. As temperature 
falls, this diffusion inhibition makes recombination more likely than annihilation with 
new partners. The Takagi pair lifetime becomes shorter, but the fluctuations resulting 
from their diffusion become biased, so their contribution to dielectric and NMR 
relaxation weakens. In this same temperature range, the Takagi group population 
temperature dependence grows weak, because many Takagi groups become stranded, 
finding no annihilation partner. These stranded groups cause relaxation in their vicinity 
but no relaxation elsewhere, so the crystal is nonergodic on the time scale concerned. 2D 
NMR shows the mean (averaged over the crystal) time required for a given nucleus to see 
a change caused by such diffusion. 1 We propose 4D NMR 2 to test whether this time is 
different in different crystal nanoregions because of nonergodicity, and to find the time 
required to restore ergodicity. Another proposed application of 4D NMR is for proton 
glass crystals with small x, for which ferroelectric nanoclusters coexist with disordered 
proton glass regions at low temperatures. 3 We expect that 4D NMR can find the size and 
lifetime of such polar nanoclusters. Such experiments could also study polar nanocluster 
size and lifetime in relaxor ferroelectric and related high-piezoelectric-strain crystals. 
! J. Dolinsek, B. Zalar, and R. Blinc, Phys. Rev. B 50, 805-821 (1994). 
2 K. Schmidt-Rohr and H.W. Spiess, Phys. Rev. Lett. 66, 3020-3023 (1991). 
3 Z. Trybula, V.H. Schmidt, and J.E. Drumheller, Phys. Rev. B 43, 1287-1289 (1991). 
♦Work supported by NSF Grant DMR-9805272. 



PYROELECTRIC SCANNING MICROSCOPY OF FERROELECTRIC POLYMER FILMS 



Stephen Ducharme (a), Brad Petersen (a), Candace Bacon (a), L. M. 
Blinov (b) , V. M . Fridkin (b,a), S. P. Palto (b) , N. N. Petukhova 
(b), S. G. Yudin (b) 

a) Department of Physics and Astronomy, Center for Materials Research 
and Analysis University of Nebraska, Lincoln, Nebraska, USA, 
sducharmel@unl . edu 

b) Institute of Crystallography, Russian Academy of Sciences, Moscow, RUSSIA 

We have developed a new technique for imaging the polarization of 
ferroelectric films and devices-Pyroelectric Scanning Microscopy 
(PSM) . PSM consists of scanning a focused laser beam across a 
capacitor made of a pyroelectric material. The laser beam is 
modulated, causing a small modulation in the local film temperature 
and generating a small ac current due to the pyroelectric effect. 
Since the pyroelectric current is proportional to the polarization, 
this technique allows us to image the polarization of the film. 
Ordinary microscopy affords a resolution of about 1 micrometer while 
near- field scanning microscopy offers resolutions of 50 nm or better. 
PSM images of ferroelectric polymer Langmuir-Blodgett films only a 
few nm thick reveal film inhomogeneities and edge-enhancement effects 
at intermediate states in the hysteresis loops. The PSM technique can 
also produce time-resolved movies of dynamic processes like switching 
by recording switching transients at each image spot and combining 
the time-space data into image frames, provided that the switching 
transients are repeatable [1] . The dynamic PS method allows us to 
monitor inhomogeneity in dynamics and should allow us to record 
domain motion. See our first Polarization Switching Movie-colorized, 
uncut, and uncensored-at 

physics .unl.edu/directory/ducharme/ducharme.html . 

Work at the University of Nebraska was supported by the USA National 
Science Foundation, the USA Office of Naval Research, the Nebraska 
Research Initiative, and the J. A. Woollam Company. Work at the 
Institute of Crystallography was supported by the Russian Foundation 
for Basic Research (#99-02-16484) and the Inco-Copernicus Programme 
(#IC15-CT96-0744) . 

1.' S. Ducharme, M . Bai, M. Poulsen, S. Adenwalla, S. P. Palto, L. M. 
Blinov, V. M. Fridkin, Ferroelectrics to be published (2001). 
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Characterising high-pressure phase transitions by X-ray Diffraction. 

R.J. Angel, N.L. Ross (Crystallography Laboratory, Dept. Geological Sciences, Virginia Tech, 
Blacksburg VA 24060, USA). 

While structural phase transitions as a function of temperature have been studied extensively, 
relatively little work has been performed at high pressures. The application of pressure to such 
systems is expected to strongly renormalise the order parameter behaviour and to suppress any 
dynamic behaviour associated either with intermediate "flip mode" phases or the high-temperature 
para phases. In addition, pressure couples linearly with the strain in the expression of the free 
energy of the solid, through the PV term. Thus, unlike temperature, pressure can drive transitions 
directly through the strain, or by inducing changes in the elastic constants. 

Recent advances in the experimental techniques for single-crystal high-pressure X-ray diffraction 
now allow lattice parameters to be determined routinely to a precision of -1 part in 40,000. The 
use of internal diffraction standards to measure pressure yields a precision in pressure 
measurement of better than 0.01 GPa over the pressure range of 0-10 GPa. These advances 
therefore open up the opportunity of characterising structural phase transitions at high pressures 
through the determination of the evolution of the spontaneous strain accompanying the transition. 
Because the coupling of the spontaneous strain to the primary order parameter of the phase 
transition is governed by symmetry rules, lattice parameter measurements can reveal the order 
parameter behaviour of the system. Examples of recent work on high-pressure phase transitions in 
titanites, melilites and lead phosphate will be presented. 



Nonexponential Response in Piezoelectric PVDF 



G. W. Bohannan 
MSU Physics 
Bozeman, MT 59717 

Polymers are known for displaying nonexponential response to both electrical and 
mechanical stresses. A numerical method for integrating fractional differential 
equations has made it possible to verify the connection between a power-law ac 
permittivity function and power-law decay currents resulting from application of a 
dc step voltage. The response shows significant dependence on the duration of the 
applied pulse. The history dependence predicted by the simulations has been 
confirmed experimentally. The computational techniques should be applicable to 
a wide range of systems exhibiting nonexponential response. 



Computer Simulations of Domain Patterns in Ferroelectrics 

Wenwu Cao and Ahluwalia 
Materials Research Laboratory 
The Pennsylvania State University, University Park, PA 16802 



We study domain pattern formation in ferroelectrics based on a 2-D time-dependent 
Ginzburg Landau theory. The electrostrictive and elastic effects are included in the form 
of nonlocal long-range interactions of the polarization field, i.e., by eliminating the strain 
variables subject to the elastic compatibility constraints. The 2-D system under study has 
a square to rectangle symmetry change during the paraelectric to ferroelectric phase 
transition, which creates four equivalent polarization states in the ferroelectric phase. We 
start the simulation from a paraelectric state that is quenched below the phase transition 
temperature. Domain pattern evolutions are obtained in absence of external electric field 
from random field initiated nucleation. In our model, the depolarization field is assumed 
to be completely screened by surface charges. For the case without defects, the final 
domain patterns are twinned with only head to tail (uncharged) domain walls. However, 
for the case with randomly distributed dipolar defects, head-to-head and tail-to-tail 
charged walls are also observed. These results are in accordance with recent experiments 
on PMN-PT and PZN-PT domain engineered single crystals, for which charged domain 
walls have been reported. 



First-Principles Effective Hamiltonians for PbSc 1 / 2 Nb 1 / 2 03 

ERIC COCKAYNE, BENJAMIN P. BURTON, 
Materials Science and Engineering Laboratory, Ceramics Division 
National Institute of Standards and Technology, Gaithersburg, MD 20899 

LAURENT BELLAICHE 
Department of Physics 
University of Arkansas, Fayetteville, AR 72701 

PbSq^Nbj^Oa (PSN) is a relaxor ferroelectric whose dielectric properties 
depend on the amount of fcc-type ordering of Sc and Nb on the perovskite 
B-site sublattice. We develop first-principles-based effective Hamiltonians 
for comparing the temperature-dependent dielectric behavior of ordered 
and disordered PSN. We compare the virtual crystal approach to modeling 
arbitrary cation configurations with the supercell approach, the latter 
based on explicit calculations for differently-ordered hypothetical PSN 
structures. We determine the sensitivities of ionic effective charges and 
interatomic forces to their environments and identify which characteristics 
are transferable between PSN structures with different cation 
configurations. 
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Multiple Aging Phenomena in Relaxor Ferroelectrics 

Eugene V, Colla, Lambert K. Chao, M. B. Weissman 
University of Illinois, Urbana, IL 61801 



The similarity of many macroscopic properties of relaxor ferroelectrics and spin 
glasses has led to consideration of the possibility of a close analogy between the 
collective frozen spin-glass state and a collective frozen state of interacting polar 
nanodomains. One key signature of spin glass states is a set of peculiar aging effects on 
the susceptibility %• We present the results of aging experiments on three familiar 
relaxors: PMN, PLZT 9/65/35 and (PMN) 0 . 7 (PT) 0 .3. These all show decreases in x' and 
X" as the sample ages over broad ranges of T. For PMN, qualitative aging results 
resembled those for spin glasses. So did the PLZT aging well below the freezing 
temperature, T g . However, PLZT aging near T g was dominated by other large effects. 
PMN-PT aging lacked qualitative memory effects found in spin glasses, resembling 
instead dirty ferromagnets. Even in PMN, the effective field scale required to erase 
memory of aging was much larger than the analogous scale in spin glasses. We have 
preliminary Barkhausen noise results on PMN, which, together with the aging, should 
help resolve the roles of collective inter-nanodomain structures and random fields in 
setting the field sensitivity of the aging. 



Electron-Nuclear Resonance Spectroscopic Evidence for Defect Pinning of 
Microdomains in Ferroelectrics : Effect of Charge and Site Substitution 

Naresh Dalai, Department of Chemistry and National High Magnetic Field 
Laboratory, Florida State University, Tallahassee, FL 32312 

and 

Dwight Viehland, Department of the Navy, Newport, RI 02841 

It has been well-established that defect-related charge imbalance and 
bonding effects play a significant role in controlling essentially all 
the properties of ferroelectric crystals and their ceramics and 
glasses. However, the details of the underlying mechanism at the 
atomic level are not yet fully understood. With the view of obtaining 
such details, we have carried out a systematic NMR and EPR study of 
this phenomenon, starting with the KH2P04 (KDP) family. We have 
introduced known quantity of structural defects, introduced by doping 
together with radiation exposure. The introduced defects have 
different sizes , as well charges, and are substituted at both the 
cation and anionic sites. Temperature dependence of spectroscopic 
splittings have been used as dynamic probes of the pinning of the 
microdomains around the defect. Electric field and isotopic effects 
on the Tc have been used to verify the proposed mechanism of the 
microdomain pinning. The presentation will focus on the experimental 
details of the defect control, dynamics measurement, as well 
theoeretical models of the observed effects. 
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SWITCHING DYNAMICS OF TWO-DIMENSIONAL FERROELECTRIC POLYMER FILMS 

! Stephen Ducharme (a), Timothy Reece (a), A. V. Sorokin (a,c), L. M. 

! Blinov (b) , V. M . Fridkin (b,a), S. P. Palto (b) , K. Verkhovskaya 

(b), N. N. Petukhova (b) , S. G. Yudin (b) 

j a) Department of Physics and Astronomy, Center for Materials Research 

' and Analysis University of Nebraska, Lincoln, Nebraska, USA, 

sducharmel@unl . edu 

b) Institute of Crystallography, Russian Academy of Sciences, Moscow, RUSSIA 

c) Department of Physics, Ivanovo State University, Ivanovo, RUSSIA 

The Ferroelectric films of vinylidene fluoride copolymers show 
j finite-size scaling of the coercive field down to thickness of 15 nm, 

I and switching at the intrinsic coercive field of about 500 MV/m in 

the 15 nm and thinner films, which are made by Langmuir-Blodgett 
deposition [1], presumably because nucleation is suppressed [2]. 
While extrinsic switching in the thicker films can be achieved in 
1/10 microsecond with exponential dependence on the switching field 
[3], intrinsic switching in the ultrathin films, 15 nm and thinner, 
[ takes 10s of seconds, with weak dependence on the switching field 

l [4]. We have succeeded in preparing films thinner than 15 nm that 

switch in microseconds with coercive fields of about 50 MV/m and 
j strong dependence of switching time on switching voltage, presumable 

j due to the introduction of new nucleation processes that are not 

suppressed by nm thickness. These films switch as quickly as 1 
microsecond at about 5 V, making the LB ferroelectric copolymer a 
J promising candidate for low-voltage non-volatile random access 

[ memories . 

Work at the University of Nebraska was supported by the USA National 
, Science Foundation, the USA Office of Naval Research, the Nebraska 

Research Initiative, and the J . A. Woollam Company. Work at the 
Institute of Crystallography was supported by the Russian Foundation 
for Basic Research (#99-02-16484) and the Inco-Copernicus Programme 
i (#IC15-CT96-0744) . 

1. S. Ducharme, V. M . Fridkin, A. Bune, S. P. Palto, L. M . Blinov, 
N. N. Petukhova, S. G. Yudin, Phys . Rev. Lett. 84, 175 (2000). 

2. A. K. Tagantsev, Integrated Ferroelectrics 16, 237 (1997). 

3. T. Furukawa, H. Matsuzaki, M . Shiina, Y. Tajitsu, Jpn. J . Appl. 
i Phys. 24, L661 (1985) . 

4. S. Ducharme, M . Bai, M. Poulsen, S. Adenwalla, S. P. Palto, L. M. 
Blinov, V. M. Fridkin, Ferroelectrics to be published (2001). 
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NANOSCALE PHENOMENA IN FERROELECTRIC THIN FILMS 



Chandan S. Ganpule 1 

Materials Research Science and Engineering Center, University of Maryland, College 

Park, MD 20742 

Ferroelectric materials are a subject of intense research as potential candidates for 
applications in non-volatile ferroelectric random access memories (FeRAM), 
piezoelectric actuators, infrared detectors, optical switches and as high dielectric constant 
materials for dynamic random access memories (DRAMs). With current trends on 
miniaturization, it becomes important that the fundamental aspects of scaling of 
ferroelectric and piezoelectric properties in these devices be studied thoroughly and its 
impact on the device reliability accessed. In this talk I will present our results on scaling 
studies of focused ion beam milled submicron ferroelectric capacitors using a variety of 
scanning probe characterization tools. 

In addition to scaling, retention of the logic state in FeRAMs needs further study. In 
an attempt to understand the nanoscale origins of back-switching of ferroelectric 
domains, we have examined the time dependent relaxation of remnant polarization in 
epitaxial lead zirconate titanate (PbZro.2Tio.gO3, PZT) ferroelectric thin films (used as a 
model system), containing a uniform 2-dimensional grid of 90° domains (c-axis in the 
plane of the film) using voltage modulated scanning force microscopy. Relaxation occurs 
through the nucleation and growth of reverse domains, which subsequently coalesce and 
consume the entire region as a function of time. Results on the effect of local curvature, 
faceting, and pinning of 180° domain walls on relaxation kinetics will be presented. This 
work is supported by National Science Foundation - Materials Research Science and 
Engineering Center (NSF-MRSEC). 



Electronic Mail: chandan@wam.umd.edu 



Ferroelectric Domains at Nanoscale 



A . Gruverman 

North Carolina State University 

Department of Materials Science and Engineering, Raleigh, NC 27695 

Investigation of nanoscale properties of ferroelectric materials, 
particularly in a thin layer form, is currently of great interest due to 
the development of a new generation of electronic devices. One of the most 
important issues in the physics of f erroelectrics is the formation and 
evolution of complex domain configurations. Direct imaging of domain 
structures in thin films and investigation of their behaviour under an 
applied electric field is essential, since it can provide a clue for 
general understanding of polarization reversal phenomena and degradation 
effects in ferroelectric films. 

Recently, it was shown that scanning force microscopy (SFM) is a 
well-suited technique both for nanoscale imaging and control of 
ferroelectric domains. The SFM imaging method is based on the detection of 
the piezoelectric vibration of the ferroelectric film excited by an 
external ac voltage applied through the metallized tip. We demonstrate 
that application of SFM can successfully address questions of both dynamic 
and static behavior of domains and domain walls. Ferroelectric domains as 
small as 20 nm in diameter can be successfully controlled in SFM. This 
technique provides unique information on basic mechanisms of polarization 
decay and fatigue in ferroelectric thin films, point-to-point recording of 
hysteresis loops at the nanometer scale, measurement of transient response 
during nanoscale ferroelectric switching as well as mapping of leakage 
current sites. 



First-Principles and semi-empirical calculations of Atoic and electronic structure for 
SrTi0 3 (100) and (110) Surfaces 

Eugene Heifets, and Ronald E. Cohen, Carnegie Institution of Washington, 5251 Broad Branch 
Rd., N.W. Washington D.C. 20015 and California Institute of Technology, MS 252-21, Pasadena 
CA 91 125 and Roberts L Eglitis, Eugene AKotomin, G.Borstel, Dept. of Physics, University of 
Osnabrueck, D-49069 Osnabrueck, Germany. 

We present and discuss results of the calculations of the surface relaxation for SrTi0 3 surfaces 
with different terminations using a wide variety of methods of computational physics, such as shell 
model as well as ab initio methods based on HF and DFT formalism. By means of semi-empirical 
shell model, the positions of atoms in 16 near-surface layers placed atop a slab of rigid ions are 
optimized. This permits us determination of a surface rumpling and surface-induced dipole 
moments (polarization) for different terminations of the (100) and (1 10) surfaces. Simultaneously, 
we performed ab initio calculations employing as Hartree-Fock , so several different DFT and 
hybrid techniques. Our shell model results for the (100) surfaces are in good agreement with both 
our ab initio calculations and LEED experiments. For the (110) SrTi0 3 surfaces O-termination is 
predicted to be the lowest in energy. 



Properties of PZT 95/5 from first principles 



Kevin Leung and Alan F. Wright, Sandia National Laboratories 
Eric Cockayne, National Institute of Standards 



The solid solution Pb Zr_0.95 Ti_0.05 (PZT 95/5) undergoes 
antiferroelectric to ferroelectric phase transition when 
temperature and pressure vary. This material is therefore 
useful in many technological applications such as power supply. 
Since large single crystal samples for the material are 
not readily available, first principles results are particularly 
valuable. We investigate various stable and metastable 
phases of both PZT 95/5 and PZ from first principles, in 
the LDA approximation, using pseudopotentials . The virtual 
crystal approximation is applied to treat the solid solution. 
The zero temperature structure for PZ obtained from first principl 
is compared to recent X-ray and neutron-scattering experiments. 
Properties such as polarization and the elastic tensor in various 
phases are computed. A modified ab initio statistical mechanical 
Hamiltonian is constructed from the potential energy surface of 
PZT 95/5, and the zero temperature structures obtained therefrom 
are compared to the ab initio structures. We also make a 
comparative study of the energetics and structures in several 
phases computed using different types of pseudopotentials. 



This work was supported by the Department of Energy under Contract 
DE-AC04-94AL85000 . Sandia is a multiprogram laboratory operated 
by Sandia Corporation, a Lockheed Martin Company, for the U.S. 
Department of Energy. 



Virtual Crystal Approximation Study of Monoclinic PZT 



Nathan R. Martin, Department of Chemistry, University of Pennsylvania 

Nicholas J. Ramer, Department of Chemistry, Long Island 
University-C. W. Post Campus 

and Andrew M. Rappe, Department of Chemistry, University of Pennsylvania 

Recently, a new monoclinic phase (space group {\em Cm}) of 
Pb(Zr$_{l-x}$Ti$_x$)0$_3$ (PZT) has been experimentally determined at 
low temperatures for compositions between $x=0.45$ and $x=0.5$ [1]. 
It has been suggested [2] that this phase may explain the large 
electromechanical response found in PZT ceramics at these same 
compositions, since its location in the phase diagram is adjacent to 
the low- temperature rhombohedral ({\em R3m} ) , high-temperature 
rhombohedral ({\em R3c}) and tetragonal ({\em P4mm} ) phases. 

In this study, we use a first-principles virtual crystal approximation 
approach [3] to study these four phases in PZT for compositions 
between $x=0.4$ and $x=0.55$. We determine the energy differences 
between these phases and construct a theoretical phase diagram at low 
temperatures which can be compared to the experimentally-determined 
diagram. 

[1] B. Noheda, D. E. Cox, G. Shirane, J. A. Gonzalo, L. E. Cross and 

S-E. Park, Appl. Phys . Lett. 74, 2059 (1999); B. Noheda, 

J. A. Gonzalo, L. E. Cross, R. Guo, S-E. Park, D. E. Cox and 

G. Shirane, Phys. Rev. B 61, 8687 (2000). 

[2] R. Guo, L. E. Cross, S-E. Park, B. Noheda, D. E. Cox and G. Shirane, 
Phys. Rev. Lett. 84, 5423 (2000). 

[3] N. J. Ramer and A. M. Rappe, Phys. Rev. B 62, 743 (2000). 



Molecular Dynamic Simulations of MgSi03 Perovskite 
and the Composition of the Lower Mantle 

Frederic C. Marton and Ronald E . Cohen 

Geophysical Laboratory and Center for High Pressure Research, 
Carnegie Institution of Washington 

Orthorhombic (Pbnm) (Mg,Fe)Si03 perovskite accounts for perhaps 
two- thirds of the volume of the Earth's lower mantle, so an 
understanding of its thermal equation of state (EoS) and thermoelastic 
properties is vital for our understanding of this region. Static 
compression experiments, however, have covered only the upper third of 
the lower mantle's geothermal gradient and elastic wave velocity 
measurements have only been carried out at pressures of up to a few GPa 
and temperatures < 1000 K. Using molecular dynamics simulations at 
pressures and temperatures of the lower mantle, we investigated the EoS 
and elasticity of MgSi03 perovskite. The simulations were done using the 
non-empirical VIB interatomic potential. This is a many-body potential 
obtained by overlapping spherical ions. Nominally charged ions give 
semi-quantitative, but less accurate results than desired for MgSi03 
perovskite. By varying the ionic charges, we found it possible to obtain 
excellent agreement with the much more computationally intensive 
first-principles LAPW method. 

Our P-V-T results were fitted to a thermal EoS of the form P(V,T) = 
P0(V,T0) + Delta Pth(T), where TO = 300 K and P0 is the isothermal 
Universal EoS [Vinet et al . , 1987]. The thermal pressure Delta Pth can 
be represented by a linear relationship Delta Pth = a + bT [Anderson, 
1980] . Using data calculated at pressures of 0-140 GPa and temperatures 
of 300-3000 K, we find V0 = 165.40 A A 3, KT0 = 273 GPa, K • TO = 3.86, a = 
-1.99 GPa, and b = 0.00664 GPa/K. Elastic constants were calculated via 
the lattice fluctuation formula of Parrinello and Rahman [1982]. 
Calculations at 300 K and 0 GPa give KT = 269 GPa and mu = 159 GPa, 
within 2% and 10%, respectively, of the experimentally derived values of 
264 and 177 GPa of Yeganeh-Haeri [1994] and in excellent agreement with 
P-V-T EoS value of 273 GPa. Results for KT at temperatures of 1000-3000 
K and pressures of 25-140 GPa agree within 0.1-5% with values obtained 
via our P-V-T EoS. 

By fixing V0 to the experimentally determined value of 162.49 A~3 [Mao 
et al., 1991] and calculating density and seismic velocity profiles 
along a lower mantle geothermal gradient, we find that the lower mantle 
cannot consist solely of (Mg,Fe)Si03 perovskite with XMg ranging from 
0.9-1.0. Using simplified pyrolitic compositions of 67 mol% perovskite 
(XMg = 0.93-0.96) and 33 mol% magnesiowusti te [(Mg,Fe)0] (XMg = 
0.82-0.86), however, we obtained density and velocity profiles that are 
in excellent agreement with seismological models. 



Growth and control of Domain Structure of epitaxial PbZro.2Tio.sO3 films grown on 

vicinal (001) SrTi0 3 

V.Nagarajan, C.S Ganpule, S.P.Alpay, A.Roytburd and R.Ramesh. 

Highly tetragonal, epitaxial PZT films with a nominal composition of 
PbZro.2Tio.8O3 exhibit a 2-dimensional grid of 90° domains (a-domains, i.e., c-axis in the 
plane of the film). Our previous studies have revealed that the 90° domains are 
preferential sites for the nucleation of 180° reverse domains during polarization switching 
and relaxation. Furthermore, we have observed that this array of 90° domains effectively 
isolate neighboring c-axis oriented regions. Therefore, we are studying approaches to 
control the spacing and periodicity of the 90° domains. Such self-assembled arrays of 
periodic domain structures can form the templates for novel memory arrays. In this 
paper, we report on the use of vicinally cut [along [100], [010], and [1 10] directions in 
the substrate plane] single crystal substrates to control the 90° domain formation. 
Epitaxial thin films have been deposited by pulsed laser deposition onto the vicinal 
substrates, with epitaxial conducting oxide bottom electrodes (LSCO and SRO). We 
have been able to control the nucleation of the 90° twins to occur preferentially at the 
steps on the substrate. We show that the orientation of these domains can be controlled 
such that they exhibit only 2 of the 4 possible variants. By using 4-circle x-ray 
diffraction, TEM and Electric Force Microscopy (EFM) we have investigated the 
structural and electrical properties of these artificially engineered structures. The control 
of such structures as a function of film thickness and substrate miscut orientation will be 
presented. 



Coexisting electronic and magnetic order in dense oxygen 

J. B. Neaton (Rutgers University), N. W. Ashcroft (Cornell University) 

Using density functional theory within the local spin-density 
approximation (LSDA) and also with gradient corrections, a new 
insulating phase of dense oxygen is found to be favored between 10 and 
65 GPa, a pressure range coinciding with that of the previously 
unidentified $\epsilon$-phase . Most strikingly, this phase is both 
antif erromagnetic and antif erroelectric . The phase has an 
orthorhombic 4-molecule unit cell with space group $Cmcm$ and is lower 
in enthalpy than both $Fmmm$ (ascribed to the $\delta$ -phase) and 
$C2mm$ (proposed for the $\zeta$-phase) . In addition, a small 
monoclinic distortion of $Cmcm$ also appears to be energetically 
favorable. The existence of antif erroelectricity may explain the 
pronounced increase in IR activity observed in the $ \epsilon$ -phase . 



Pseudopotential Tail Norm Conservation for Ferroelectric Modeling 
Nicholas J. Ramer, 

Department of Chemistry, Long Island University-C . W. Post Campus 
Ilya Grinberg and Andrew M. Rappe, 

Department of Chemistry, University of Pennsylvania 

Ferroelectric phenomena result from the interplay of atomic size and 
atomic polarizability . Therefore, it is vital to use pseudopotentials 
(PSPs) which accurately represent these atomic properties, to insure 
the accuracy of pseudopotential (PSP) density functional calculations 
of ferroelectric materials. The preservation of atomic size for the 
reference configuration has been formulated into our implementation of 
the virtual crystal approximation [1] and relativistic PSP [2] 
constructions. However these methods do not insure the preservation 
of the atomic size as a function of electronic configuration. 
Therefore we have recently imposed an • additional atomic criterion in 
our PSP constructions. We quantify the atomic size by examining 
$\frac{dN_i} {df_j }$ (the change in wave function tail norm $N$ for an 
atomic orbital $i$ as a function of occupation for an atomic orbital 
$j$) . We will show that the preservation of this quantity 
dramatically affects the accuracy of the resulting structural and 
elastic properties of materials in solid-state calculations. 

[1] N. J. Ramer and A. M. Rappe, Phys . Rev. B 62, 743 (2000). 



[2] I. Grinberg, N. J. Ramer, and A. M. Rappe, Phys. Rev. B 61, 1 (200 



Ferroelectric and piezoelectric properties in the presence of 
compositionally broken inversion symmetry 

Na Sai, Bernd Meyer and David Vanderbilt 
Department of Physics and Astronomy 
Rutgers University 
Piscataway, NJ 08855-0849 
nsai@physics.rutgers.edu 

We have investigated novel ferroelectric perovskite systems in which compo- 
sitional inversion symmetry is broken [1]. Our method involves the calculation 
of the density-functional total energy and Berry-phase polarization. Here we ex- 
tend the previous work by (i) focusing on the piezoelectric response of the triple- 
cell heterovalent system Ba(Ti-<5,Ti,Ti-f 5)03, and (ii) studying the strength of 
the symmetry breaking in the double-cell system (Ba,Sr)(Ti-<5,Ti+<$)03 with 
simultaneous A-site and B-site substitutions. We observe an enhanced piezo- 
electric response coefficient e33 when increasing the compositional parameter 6 
in the triple-cell system. This enhancement is quite drastic for the metastable 
minimum, but only modest for the stable minimum, corresponding respectively 
to the minority and majority wells of the ferroelectric double-well structure. In 
the double-cell system, we find that the increase in strength of symmetry break- 
ing with the chemical or concentration perturbation is dominated by the term 
linear in 6 y in contrast with the S 3 dependence found in the triple-cell system. 
A symmetry-based justification of the dominance of the linear term is provided. 



[1] N. Sai, B. Meyer and D. Vanderbilt, Phys. Rev. Lett. 84, 5636 (2000). 



Kinetic Monte Carlo Simulations of Crystal Growth in Ferroelectric Material 

M. Suewattana, C. H. Tahan, P. Larsen, Shiwei Zhang, and H. Krakauer 
Dept. of Physics, College of William and Mary 
Williamsburg, VA 23187-8795 



We study the growth process of ferroelectric materials by kinetic Monte 
Carlo simulations. An electrostatic model with long-range Coulomb interac- 
tions [1] is used to model relaxor single crystals. The growth is character- . 
ized by thermodynamic processes involving adsorption and evaporation, with 
solid-on-solid restrictions. A new algorithm is developed in order to simulate 
growth under such a model, for which existing kinetic Monte Carlo algo- 
rithms are inadequate. We study the growth rates and the order structure 
of the grown crystals as a function of temperature, chemical composition,, 
and growth orientation. Tests of our algorithm on NaCl give good results. 
We discuss our on-going simulations of the growth processes in PMN-type 
materials. 

* Supported by the Office of Naval Research. 

[1] L. Bellaiche and D. Vanderbilt, Phys. Rev. Lett. 81, 1318 (1998); 
M. Wu and H. Krakauer, http://xxx.lanl.gov/abs/cond-mat/0011063, to be 
published. 
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Pressure Dependence of the Phase Diagram for PZN/PT and PMN/PT Single Crystals 

E.L. Venturing G. A. Samara and B. A. Turtle, Sandia National Laboratories, 
V.H. Schmidt, Montana State University, and S.-E. Park, Pennsylvania State University 

We employ hydrostatic pressure to elucidate the nature of the high-temperature cubic-tetragonal 
ferroelectric (FE) transition and the lower temperature tetragonal-rhombohedral (T-R) structural transition 
in single crystals of PZN/9.5%PT and PMN/33%PT, compositions near the morphological phase boundary 
(MPB). This information provides new insight into the physics of these technologically important mixed- 
oxide materials. The T-R transition is detected in both dielectric constant and dissipation data for both 
crystals and exhibits some sample-history-dependent effects. The transition temperature shows a gradual 
increase with pressure in PZN/PT and disappears at pressures above 10 kbar where the T-R phase 
boundary meets the FE transition. For PMN/PT the T-R transition temperature is nearly independent of 
pressure up to 7 kbar. 

The FE transition temperature in PZN/PT decreases at -4.7 K/kbar from -455 K at ambient pressure; 
PMN/PT exhibits a larger pressure dependence of -5.3 K/kbar from -420 K at ambient. More importantly, 
the nature of the FE transition in PZN/PT changes to a diffuse or relaxor characteristic at pressures above 5 
kbar as shown by an increasing frequency dispersion in the dielectric constant maximum temperature. The 
change in the dynamics of the relaxor transition with pressure can be analyzed by fitting the dielectric 
constant maximum temperature versus frequency to the Vogel-Fulcher equation. The fitting parameters are 
an energy barrier and "interaction" temperature: the energy barrier increases with increasing pressure while 
the "interaction" temperature decreases -20% faster with pressure than the dielectric maximum 
temperature at a fixed frequency. These changes are consistent with an increasing isolation of shrinking 
nanopolar FE domains at higher pressures. In contrast to PZN/PT, no crossover to relaxor behavior is 
observed in the PMN/PT crystal at pressures up to 20 kbar. 

Acknowledgement: Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed 
Martin Company, for the U.S. Department of Energy under Contract No. DE-AC04-94AL85000. 



NEW POLARONIC-TYPE EXCITONS IN FERROELECTRIC OXIDES: NATURE 
AND EXPERIMENTAL MANIFESTATION 



V. S. Vikhnin. R. I. Egiitis* , E.A. Kotomin + * , and S. Kapphan + 

A. F. Ioffe Physical Technical Institute, 194021, Saint- Petersburg, Russia 
+ University of Osnabrueck, Fachbereich Physik, 49069 Osnabrueck 
* University of Riga, Kengaraga 8, 1234 Riga, Latvia 

Charge Transfer Vibronic Exiton (CTVE) in ferroelectric oxides with mixed ionic-covalent 
type of chemical bonding are correlated pairs or/and triads of electronic and hole polarons 
[1,2]. Ferroelectric oxides are good examples of such ionic-covalent bonding cases with rather 
strong charge transfer and not small vibronic interaction. Ferroelectric oxides KTa0 3 (KTO), 
KNb0 3 (KNO), and solid solutions KTN and SBN - bulk as well as thin film materials will be 
considered as examples of new type excitons of polaronic type (CTVE). Here a review of 
recent theoretical and experimental investigations of CTVE is presented. 
The important feature of CTVE is the formation of charge transfer which is self-consistent 
with accompanied lattice distortion. The computations of electronic and distorted lattiqe 
structure of CTVE in KTO-crystal was performed on the basis of semi-empirical Hartree-Fock 
calculations using the intermediate neglect of the differential overlap (INDO) method. We have 
also used the periodic large unit cell method and have extended the primitive KTO unit cell 
3x3x3 i. e. 27 times. As a result, the following structure of CTVE was obtained. The CTVE 
here has a triad structure [2] with rather strong vibronic energy lowering (~ 2,71 eV). This 
energy lowering is connected with the appearance of a hole polaron on one oxygen ion which 
has the displacement 5,2 % towards to the first active Ta-ion where the electronic polaron is 
partly located. This first Ta-ion has a displacement of 3,1 % towards the active oxygen ion. But 
the second active Ta-ion, where the electronic polaron is also partly located, has a displacement 
of 4,5 % with repulsion from the active oxygen ion in contrast to the first one. The final charge 
distribution on these three Ta-O-Ta ions is follow-ing: the first Ta-ion has + 1,74 total charge, 
while the second Ta-ion has + 2,12 total charge, and active oxygen ion has - 0,24 total charge. 
The energy level of the electronic polaron on the two Ta-ions is about ~ 0,8 eV lower than the 
bottom of the conduction band, and the ener-gy level of the hole polaron on the oxygen ion is 
about - 1 eV higher than the top of valence band of the KTO-crystal. Similar triad structure of 
CTVE was obtained in KNO [3]. 

It should be underlined that the energy of a possible electron-hole polaron recombination with 
CTVE collapse corresponds to a luminescence quanta energy equal to about ~ 2.14 eV. This 
value is in a good agreement with experimental value of the emission quanta energy of detected 
in KTO so-called «red luminescence» signal [4]. Moreover, as it was shown in [4], the «red 
luminescence» signal connects with intrinsic defects and is disappearing under the action of 
any type of heat treatment. Both these properties are in good agreement with the assumption 
that «red luminescence)) in KTO is the CTVE recombination luminescence. Note also that 
CTVE trapped by Cr or Ce impurities in SBN (S. Kapphan, Ming Gao, V.S. Vikhnin et al, 
2000, and S.Avanesyan, H.Liu, V.S. Vikhnin, 2000), or CTVE trapped by Ca-impurity in Ba,. 
x Ca x Ti0 3 (S. Kapphan, M. Wierschem, V.S. Vikhnin et al, 2000) can be responsible for 
uncommon visible range light absorption peak. Supported in part by DAAD. 
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Charge Transfer Electrostatic Model of Compositional Order in Perovskite Alloys* 

Michael Wu and Henry Krakauer 
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We use an electrostatic model including charge transfer to study com- 
positional ordering in perovskite alloys with A-site and B-site doping. The 
method [1] generalizes the approach of Bellaiche and Vanderbilt [2], who used 
nominal ionic charges. The large covalency of Pb 2+ compared to Ba 2+ is mod- 
eled by an environment dependent effective A-site charge. Monte Carlo sim- 
ulations of this model reproduce the long range compositional order of both 
Pb-based and Ba-based A(BB'B")0 3 perovskite alloys as a function of tem- 
perature and B 4+ concentration. In systems with both A- and B-site doping, 
such as (Na 1/2 La 1/2 )(Mg 1/3 Ta2/3)03 ; (Ba 1 _ I La x )(Mg (1+x)/3 Nb { 2- x )/3)03 and 
(Pbi- z La z )(Mg(i +ar )/ 3 Ta(2- x )/3)03 [3], purely electrostatic interactions repro- 
duce general trends, and charge transfer effects indicate that local structural 
relaxations can tip the balance between different B-site orderings in Pb based 
materials. 

* Supported by the Office of Naval Research. 

[1] M. Wu and H. Krakauer, http://xxx.lanl.gov/abs/cond-mat/0011063, to 
be published. 

[2] L. Bellaiche and D. Vanderbilt, Phys. Rev. Lett. 81. 1318 (1998). 
[3] L. Dupont, L. Chai ; P.K. Davis, Sol. St. Chem. Inorg. Mater. II 547. 93- 
98 (1999); J.K. Montgomery, M.A. Akbas and P.K. Davies. J. Am. Ceram. 
Soc. 82 : 3481 (1999). 
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